THE DENSITY OF THE ISE 
AND LOCAL LIMIT LAWS FOR EMBEDDED TREES 

MIREILLE BOUSQUET-MELOU AND SVANTE JANSON 

Abstract. It has been known for a few years that the occupation measure of 
several models of embedded trees converges, after a suitable normalization, to 
the random measure called ISE (Integrated SuperBrownian Excursion). Here, we 
prove a local version of this result: ISE has a (random) Holder continuous density, 
and the vertical profile of embedded trees converges to this density, at least for 
some such trees. 

As a consequence, we derive a formula for the distribution of the density of 
ISE at a given point. This follows from earlier results by Bousquet-Melou on 
convergence of the vertical profile at a fixed point. 

We also provide a recurrence relation defining the moments of the (random) 
moments of ISE. 



1. Introduction 

We consider some families of random labelled trees; the labels will be integers 
(positive or negative). Our main case is binary trees, where each node is labelled 
with the difference between the number of right steps and the number of left steps 
occurring in the path from the root to the node. In particular, the root has label 0, 
and the labels of two adjacent nodes differ by dbl. Note that the label of each node 
is simply its abscissa, if we draw the tree in the plane in such a way that the right 
[left] child of a node lies one unit to the right [left] of its parent. We call this the 
natural labelling of a binary tree. 

Given a labelled tree T, let X(j; T) be the number of nodes in T with label j; the 
sequence (X(j; T))|2 = _ 00 is the vertical profile of the tree (Figure^). 

Let T n be a random binary tree with n nodes with the uniform distribution, and 
let X n (j) := X(j;T n ) be its vertical profile. It was shown by Marckert |23l Theorem 
5] that the (random) distribution of the labels in the tree converges, after appropri- 
ate normalization, to the ISE (integrated superbrownian excursion) introduced by 
Aldous [I], see also [7j. The ISE is a random probability measure; to emphasize this 
we will usually write it as /x ISE . (Actually, the result in [2S] is stated for complete 
binary trees, i.e. binary trees where each node has either or 2 children, but the 
result transfers immediately by considering internal nodes only; see at the end of 
Section |H] for details.) Marckert's result can be stated as follows, where 7 := 2 -1 / 4 , 
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Figure 1. A binary tree having vertical profile [2, 2, 4; 2, 1, 1]. 

£{v) denotes the label of v, and 5 X is the Dirac measure at x, 
1 °° 1 d 

- X *(j) S Tn-V*j = ~Y^ S Tn-V*l(v) > Mise, n -»■ 00, (1.1) 

j=-oo v£T n 

with convergence in the space of probability measures on R. For complete binary 
trees, the result is the same, except that now 7=1 (and n has to be odd). 

Our first main result is a local version of (|1.1[) . showing that the vertical profile 
of random binary trees, properly normalized, converges to the density / ISE of /Xise! 
see Section 01 for details. Our second result consists of a recurrence relation that 
characterizes the joint law of the moments of the ISE. 

Remark 1.1. Different normalizations of /x ISE are used in the literature. We use the 
normalization of [I], also used by e.g. (2H]. The normalization in [3(25 differs by a 
scale factor 2 1 / 4 . 

Our local limit result actually holds for other families of labelled (or: embedded) 
trees too. Indeed, the random measure /U ISE arises naturally as a limit for embedded 
trees in the following way [I]. Let T n be a random conditioned Galton-Watson tree 
with n nodes, i.e. a random tree obtained as the family tree of a Galton-Watson 
process conditioned on a given total population of n. (See e.g. [21 El for details, 
and recall that this includes e.g. binary trees, complete binary trees, plane trees and 
labelled (=Cayley) trees. These random trees are also known as simply generated 
trees.) The Galton-Watson process is defined using an offspring distribution; let 
£ denote a random variable with this distribution. We assume, as usual, E£ = 1 
(the Galton-Watson process is critical) and < <r| := Var£ < 00. Assign i.i.d. 
random variables rj e to the edges of T n . We regard r? e as the displacement from 
one endpoint of the edge e to the other, in the direction from parent to child; this 
gives a labelling of the nodes such that the root has label and each other node 
v has label t(v) := £(v') + r] vv >, where v' is the parent of v. For the purposes of 
this paper, we assume rj e to be integer valued. We further assume Er/ e = and 
< <rjj := Var7/ e < 00. Define X and X n as was done above for binary trees with 

their natural labelling. Then 1)1.1(1 holds, with 7 := o~ x o^ 2 0], see also |17j . 

We conjecture that a local version of (|1.1[) holds in this generality, provided r/ e is 
not supported on a subgroup dTL of the integers with d > 2, but we will only prove 
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this for two special cases, viz. random plane trees with r\ e uniformly distributed on 
either {±1} or { — 1,0, 1}, see Theorem 13,61 

We state in Section [2] some properties of the (random) density function / ISB , in 
particular that it exists. The proofs are given in Section El after some preliminaries on 
the Brownian snake and the Brownian CRT (continuum random tree) in Sectional 
Our results on the local limit law are stated in Section Inland proved in Sections IT hIOI 
Some further computations of (mixed) moments of the density /ise(A) are given in 
Section ^2 Our results on moments of /i ISE are stated in Section 0] and proved in 
Section [121 

All unspecified limits below are as n — > oo. We will use C and c with various 
subscripts to denote various positive constants, not depending on n or other variables; 
for constants depending on a parameter we use C{a) and so on. 

Acknowledgement. This research was mainly done at the Mittag-Leffler Institute, 
Djursholm, Sweden, during the semester on Algebraic Combinatorics. MBM was 
partially supported by the European Commission's IHRP Programme, grant HPRN- 
CT-2001-00272, "Algebraic Combinatorics in Europe". 

We further thank Ingemar Kaj and Jean- Francois Le Gall for helpful comments. 

2. The density of the ISE 

It is no surprise that the random measure fi lSE is absolutely continuous; the fol- 
lowing theorem may well be known to experts, but we have not found an explicit 
reference. (Related results for super-Brownian motion have been given by |19 |l26[ l3"0]. 
It seems to be possible but non-trivial to derive the existence of a density for ISE 
from these results.) We give a proof in Sectional 

Theorem 2.1. ISE has a Holder continuous density. In other words, there exists a 
continuous stochastic process /ise(^); — oo < x < oo, such that dfi lSE (x) = f lSE (x)dx. 
Moreover, the random function fis E (x) has a.s. the following properties: 

(i) /ise has compact support: sup{|x| : fi SE (x) / 0} < oo; 

(ii) / ISE is Holder (a) -continuous for every a < 1; 

(iii) /ise has a derivative f( SE (x) a.e. and in distribution sense, and // SE € LP(dx) 
for every p with 2 < p < oo . 

Of course, the support of / ISE is random; (i) says that there exists a random 
M < oo such that fis E (x) = for \x\ > M, but no deterministic M will do. 

Remark 2.2. More precisely, the proof in Section El shows that / ISE belongs to the 
generalized Sobolev space L 2,a for any a < 3/2. Loosely speaking, / ISB thus has 
"a derivatives in L 2 " for every real a < 3/2. 

Parts (ii) and (iii) of Theorem 12 . 1 1 come close to showing that / ISE has a continuous 
derivative, but we have not been able to prove it. Indeed, it seems likely that the 
(fractional) derivatives in L 2 asserted by Remark 12 .21 are continuous. Hence we make 
the following conjecture. 

Conjecture 2.3. The density f lSE has a.s. a continuous derivative, but not a second 
derivative. 
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The marginal distributions of / ISE , i.e. the distributions of /i S e(A) for fixed A, will 
be described in Corollaries 13.31 and 13.41 Moments and mixed moments of /ise(A) will 
be computed in Section fTTl 



3. Local limit results 

Our main result is the following local limit result for naturally embedded random 
binary trees, conjectured in Q. 

We let X n {x) denote the function obtained by extending X n (j) to arbitrary real 
arguments by linear interpolation; thus X n (j) = X n {j) for every integer j, and X n 
is linear on each interval [j, j + 1] . 

Co(K) denotes, as usual, the Banach space of continuous functions on R that 
tend to at ±oo. We equip Co(M) with the usual uniform topology defined by the 
supremum norm. 

Recall that we have defined the constant 7 as 2~ 1 / 4 for binary trees and 1 for 
complete binary trees. 

Theorem 3.1. Consider random binary trees or random complete binary trees with 
their natural labelling. Then, as n — > 00, 

VV/ 4 X n ( 7 -V/ 4 *) / ISE (x), (3.1) 

in the space Cq(M) with the usual uniform topology. Equivalently, 

n-^X^n^x) -±> 7 /.se(7^)- (3-2) 

Note that the functions on the left-hand sides of (|3.1|) and (|3.2|) are density func- 
tions, i.e. non-negative functions with integral 1. Proofs will be given in Sections 

□EH 

Corollary 3.2. For random binary trees or random complete binary trees with 
their natural labelling, if n — > 00 and jn/n 1 ^ — > x, where —00 < x < 00, then 



n 



- 3/4 x n (j n )^ 7 f 1SE ( r x). 



It follows by combining this with results in Bousquet-Melou [7j that the marginal 
distributions of / ISE are as conjectured there. 

Corollary 3.3. For every real x, the distribution o// ISE (x) is given by the moment 
generating function 

Ee a/ ISE (x) = L(2 -V4| x | )2 -l/4 a)) | a | < 2 2+l/4 3 -l/2 j 

where 

t, ^ 1 48 f A(a/v 3 )e- 2xv 5v a, 
L(x, a) :=! + —= \ 1 v 5 e v dv, x > 0, 

W 71 Jr \i- + A{a/v A )e j 

A(y) = A is the unique solution of 

y (1 + Af 
24 I- A 
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satisfying A(0) = 0, and the integral is taken over 

T = {1 - te~ in/4 ,t £ (-oo, 0]} U {1 + te m/A , t £ [0, oo)}. 

In particular, the density at x = has a simple law. (See again 0.) 

Corollary 3.4. /i SE (0) has the same distribution as 2 1 / 4 3 -1 T~ 1 / 2 , where T is a 
positive 2/3-stable variable with Laplace transform Ee _iT = e - * 2 3 . 
Hence /ise(0) has the moments 

Ef 1SE (o) r = y^ r^+ir ~ 4/3 <r< °°- 

As said in the introduction, we conjecture that the local limit results hold also for 
conditioned Galton-Watson trees with random labellings defined by i.i.d. random 
increments r\ e along the edges; a precise formulation is as follows. Recall that the 
span of rj e is the largest integer d > 1 such that rj e a.s. is a multiple of d. 

Conjecture 3.5. Consider a random conditioned Galton-Watson tree T n with a 
random labelling defined as above by integer valued random variables n e with mean 
0, finite variance > and span 1. Then, the conclusions (|3,lj) and (j3.2j) of 

Theorem 13.11 hold, with 7 := o-~ l a 1 ^ 2 . 

If this conjecture holds, the conclusion of Corollary 13.21 holds too. 
As said in the introduction, we can prove the conjecture in two special cases, both 
considered in [7j. 

Theorem 3.6. Conjecture 13.51 holds if T n is a random plane tree and r] e is uniformly 
chosen at random from {±1} or from { — 1,0,1}. 

For these two cases, <r| = 2 and Var?? e = 1 and 2/3; hence 7 = 2 1 / 4 and 7 = 
2~ 1 / 4 3 1 / 2 , respectively. 

Remark 3.7. It follows from the proof of Theorem 13.11 in Section [7| that to prove 
Conjecture 13.51 in further cases, it suffices to prove the estimate in Lemma 17.31 

4. The moments of ISE 

Let T n be a random binary tree with n nodes, and let \i n be the following (random) 
probability distribution: 

^ n = n ^2 5 (2n)-V4%)- C 4 - 1 ) 
veT n 

As recalled in the introduction, [i n converges to fi mE . The ith moment of /i n , denoted 
mi tn , is itself a random variable: 

m i<n = -J2 ^ny'/HivY = 2" i / 4 n~ 1 - i / 4 ^ ( ( V Y- 
n v£T n veT n 

We shall prove that the sequence mi,n, ffJ2,n, • ■ ■ converges in distribution to the 
sequence mi, m 2 , ... of moments of ISE, and compute the joint moments of the m^: 

E(mfmf 2 ---m^), 
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for all (fixed) values of px,P2, ■ ■ ■ ,Pr- The moments of the mj, being the moments of 
the moments (of fii SE ) should probably be called the grand-moments of fii SE . Note 
that the grand-moments of a random probability measure, provided they do not grow 
too quickly, determine the distribution of the sequence of moments of the measure, 
and thus the distribution of the random measure. 

In order to state our result, we introduce some notation. A partition A of an 
integer k is a sequence (Ai, . . . , X p ) of non-decreasing positive integers summing to 
k. The value k is called the weight of A, also denoted k = |A|. For instance, 
A = (1,1,3,4) is a partition of k = 9. The A« are called the parts of A. We 
shall also use extended partitions, in which the positivity condition on the parts is 
relaxed by simply requiring that Aj is non-negative. Hence A = (0,0, 1, 1,3,4) is an 
extended partition of 9. The union o~L)t of two extended partitions a = (a\, . . . , o~ p ) 
and r = (n, . . . , T q ) is obtained by reordering the sequence (o"i, . . . , cr p ,n, . . . , T q ). 
For any p-tuple (a\, . . . , a p ) of non-negative integers, we denote by a the extended 
partition obtained by reordering the Uj. Given two p-tuples a and A, we write a < A 
if o"j < Aj for 1 < i < p. 

We shall denote 

m in m 2 n ' ' ' m r,n := m \,n and rn^m^ 2 ■ ■ ■ mP T r := m\ (4-2) 

where A = 1 P1 2 P2 • • • is the partition having p\ parts equal to 1, p2 parts equal to 2 
and so on. The value of E(m^) will be expressed in terms of a rational number c\, 
which we actually define for any extended partition A = (Ai, . . . , A p ). The definition 
works by induction on p + |A| as follows: 

• C0 = -2, 

• c\ = if |A| is odd, 

• c A = (p + |A|/4 - 3/2)c v if Ai = 0, with A' = (A 2 , . . . , Ap), 

• if Ai > 0, 




where J = [p]\I, A/ = (A^, X i2 , . . . , A ir ) if J = {ii, . . . , i r } with 1 < h < ■ ■ ■ < 
i r < p, the second sum runs over all non-negative p-tuples a (not necessarily 
partitions) satisfying the two required conditions, and ( cr ) = nf=i 

Theorem 4.1. As n — > oo, the moments in\ n , rfi2,m ■ ■ ■ "/ the occupation measure of 
binary trees converge jointly in distribution to the moments m\,m2, ■ ■ ■ of ISE. The 
convergence of moments holds as well, and for all partitions A, the joint \-moment 
of the random variables rrii tn , defined by (|4.2|) . satisfies 

^( m x,n) = E(wi^) =0 if |A| is odd, 

and otherwise 

E ( ">„,„) ^ E <^> - nZm-m (44) 

where the number c\ is defined just above. 
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The vanishing of E(m> jn ) when |A| is odd, is a straightforward consequence of the 
symmetry of T n . The proof is given in Section IT2l 

Example 4.2 (The average moments of ISE). When A has a single part, equal 
to 2k for k > 1, the above recurrence relation gives c/2k) = k{2k — l)c/2fc_2), together 
with the initial condition cm\ = (1 — 3/2)c0 = 1. Hence the mean of the 2A:th moment 
of the random probability \x n satisfies 

E(m 2k n ) E(m 2k ) = J^ffllZgj , (4.5) 
V lk,n) \ 2k J 2 3fc / 2 r((fc + l)/2) 

Example 4.3 (The moments of the average of ISE). Let ?7ii n denote the mean 
of \x n . Then E(m^ n ) = E(mf) = if p is odd, while 

,2fc\ . th./ 2fc\ a k Y{l/2) 



E (mf n ) E (ml 1 



2 k / 2 T((5k-l)/2) , 



where ao = —2 and for > 0, 
fe-i 



i=i ^ 

Indeed, = ca, where A = l 2k , and the recurrence relation 1)4.3)1 translates into 
the above recursive definition of a^. Note that each 2fe-tuple a occurring in 1)4.3)1 
contains 2 coefficients equal to zero, so that we also use the part of the definition 
of c\ that deals with the case Ai = 0. This value of E(m 2k ) was already obtained 
in E 



Example 4.4 (The first two moments of ISE). Let us finally work out the joint 
distribution of m\ and m 2 . We have E{m 2k ^ l m 2 n ) = E(m| fe+1 m2) = 0, and 

EK^mU - E(mf mf) - ^ r(1/2) 



2(*+*)/2r((5fc + 3£- l)/2)' 
where ao = — 2 and the t are determined by induction on k + £: 



(0,0)<(i,i)<(fe,Q 



Qi) (j) ai d a k-i/-j + 2 ^ - l)afe+i/-2 



+ -fc(2jfe - l)(5Jfe + 31 - 4)(5fc + 3£ - 6)a fc _ M + -(4fc + l)£(5k + 3£ - 4)a k/ -i. 

Here, = c\ with A = l 2k 2 e . In the right-hand side of the equation, the second 
(resp. third, fourth) term corresponds to the case a = \ 2k + 2 2^~ 2 (resp. a = 2 l 2k ~ 2 2^, 
a = 01 2fc 2^ _1 ). The last case occurs both when we replace a part of A equal to 2 by 
a zero part, and when we decrease by 1 a part equal to 1 and a part equal to 2. Of 
course, this generalizes Example 14.31 f which corresponds to t = 0). 

It seems likely that Theorem l4. ll extends to randomly labelled conditioned Galton- 
Watson trees as in Conjecture 13.5) at least under some moment conditions on £ and 
r) e , where the measure fj, n is defined by the left-hand side of 1)1.1)1 and, as usual, 
7 := a~ l a l J 2 . We show this for the special case in Theorem 13.61 
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Theorem 4.5. If T n is a random plane tree and t] e is uniformly chosen at random 
from {±1} or from { — 1,0,1}, then the conclusions of Theorem 14.11 hold, where the 
measure [i n is defined by . with 7 = 2 1 / 4 and 7 = 2~ 1 / 4 3 1 / 2 , respectively. 

To conclude this section, we want to underline briefly a similarity between the 
density of ISE and the local time of the (normalized) Brownian excursion. In fact, 
Theorem 12 . 1 1 shows that the vertical profile of a random binary tree converges, after 
suitable rescaling, to the density of the ISE. Similarly, as shown by Drmota and 
Gittenberger |10| . the horizontal profile converges to the local time of the Brownian 
excursion. 

We can develop this analogy to grand- moments as follows. Consider again the 
random binary tree T n . Let d(v) denote the depth (the distance from the root) of 
vertex v and define the probability measure 



- ^2 ^-Wn-i/Mv) > ( 4 - 6 ) 



l 2- 3 /2n-VM„) , 

veT„ 

describing the horizontal profile, i.e. the distribution of the depths (after rescaling). 

It is known, as an immediate consequence of Aldous [UO], that as n — ► 00, u n 
^excj where the random probability measure u EXC is the occupation measure of the 
Brownian excursion, and thus has the local time of the Brownian excursion as density. 

The similarity with the vertical profile and ISE is obvious, and we adopt below 
the same notation as before (Trii in ,m\ n , etc.) for the moments of v n and i^xc- I n 
particular, now rrij := f x l du EXC (x) = Jq e(t) 1 dt, where e(t) is a Brownian excursion. 
Then a result similar to Theorem 14. II holds: 

Theorem 4.6. As n — ► 00, the moments mi fi ,fn2 in ,... of the horizontal profile 
(depth distribution) measure v n converge jointly in distribution to the moments 
mi,m2,- ■■ of z^exc- The convergence of moments holds as well, and for all parti- 
tions \, the joint X-moment of the random variables I7ii, n , defined by (|4.2(1 . satisfies 

9— 3|A|/2 j tYI /9\ 

E(m A n ) E(m A ) = 1 ,\ , (4.7) 

V \nj \ \) r(p+ |A|/2 — 1/2)' V ; 

where the number d\ is defined by 

• d$ = -2, 

• d x = (p + |A|/2 - 3/2)d A ' ifM = 0, with X' = (A 2 ,...,A P ), 

• ifXx> 0, 

dx = \ E d ^ d \i+ J2 (TK> ( 4 - 8 ) 

m^iQp] «t<a,|<t|=|a|-i v j 

with the same notation as in (|4.3|) . 

The proof is very similar to the proof of Theorem 14. H but simpler. Note that the 
binomial coefficient (*) is simply equal to one of the Aj. The proof is sketched at 
the end of Section IT2l 

The grand-moments in ()4.7|) have been computed by a different method by Richard 
|2*7) : a special case (moments of mi and 7712) is given by Nguyen The [2*H] . and the 
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moments of the Brownian excursion area mi were found already by Louchard |22| . 
see also and The grand- moments in ([4.7ft . as well as the grand-moments in 
(|4.4j) above, can also be derived by the method of ^3 Section 5], which is related to 
the method used here but phrased in different terms. (Presumably, the method of 
|27] too applies to (JOj) as well.) 

Again, the same result holds for random plane trees as well, provided we change 
the scale factor 2~ 3 / 2 in to 2 -1 / 2 . 

Remark 4.7. A Dyck path of length 2n is a 1-dimensional walk starting and ending 
at 0, taking steps in {— 1,+1}, and never reaching a negative position. There is a 
well-known correspondence between plane trees with n + 1 vertices and Dyck paths 
of length 2n, where the Dyck path gives the depths of the vertices along the depth- 
first walk on the tree. It follows easily that Theorem 14.61 holds for moments of a 
uniformly chosen random Dyck path w n of length 2n too defined by 

- 2n 

mk > n := ^5Z( 2n )~ fe/2 ^( i ) fc ; 

i=l 

this has previously been shown by Richard |27j . 

Remark 4.8. It is possible to use our methods to obtain results on grand-moments 
of the vertical and horizontal occupation measures together, and thus on the joint 
distribution of the vertical and horizontal profile, and also on the asymptotic distri- 
bution of the pair of labels (£(v), d(v)) . We leave this to the reader. 

5. The Brownian snake and CRT 

5.1. The Brownian snake. We begin by recalling the definition of the Brownian 
snake, see Le Gall j201 Chapter IV] or Le Gall and Weill [2^ f° r further details; see 
also |161 Section 4.1]. Let the lifetime, be £ := 2B CK , where B ex is a Brownian 
excursion on [0, 1]. (In general, the lifetime £ might be any (locally) Holder continu- 
ous non-negative stochastic process on some interval /; in other contexts, £ is often 
taken to be reflected Brownian motion on [0, oo) [20]-) Let, for s,t € [0, 1], 

m(s,t;£) := min{((u) : u £ [s,t]} when s < t, 

and m(s, t; £) := m(t, s; £) when s > t. The Brownian snake with lifetime £ then can 
be defined as the continuous stochastic process W(s, t) on [0, 1] x [0, oo) such that, 
conditioned on £, W is Gaussian with mean and covariances 

Cov (W(s 1 ,t 1 )W(s 2 ,t 2 ) | C) =min(ti,t 2 ,m(si,s 2 ;C))- 

We have defined the Brownian snake as a random field with two parameters, but we 
are really only interested in the specialization W(s) := W(s,£(s)), s E [0,1]; this 
stochastic process is called the head of the Brownian snake. (In fact, it is easily 
seen that the pair (£, W) determines W; see further [21] •) Conditioned on £, W is a 
Gaussian process on [0, 1] with mean and covariances K(W(s)W(t) \ Q = m(s, t; £)• 
Consequently, still conditioned on £, W{s) — W(t) has a normal distribution with 
mean and variance 

Var(VF(s) - W(t) | C) = a 2 (s, t; () := ((s) + £(t) - 2m(s, t; ()■ (5.1) 
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The random probability measure /i ISB can be defined as the occupation measure of 
the process W, see j2Q], |2l] and the next subsection. Hence, / ISB is the occupation 
density of W, also called its local time. 

5.2. Brownian CRT. The Brownian CRT (continuum random tree) was intro- 
duced by Aldous ^ [H [3] as a natural limit of rescaled finite random trees. It is a 
random compact metric space that is a topological tree in the sense that every pair 
of points x,y are connected by a unique path (homeomorphic to [0,1]), and that 
path has length d(x,y). We let here and later d denote the metric. The Brownian 
CRT is further equipped with a probability measure v, which gives a meaning to "a 
random node" in the CRT. 

One of Aldous's characterizations of the Brownian CRT PI [3J Q] uses the distri- 
bution of the shape and edge lengths of the spanning subtree Rk spanned by the 
root o and k independent random nodes X\, . . . , in the tree. (Here k is an arbi- 
trary positive integer.) Then a.s., the subtree Rk admits the root and X\, . . . , X^, as 
leaves, and has exactly k — 1 internal nodes, all of of degree 3; the leaves are labelled 
but not the internal nodes. If we ignore the edge lengths (which are positive real 
numbers), there are (2k — 3)!! possible "shapes" of Rk', for each shape we number 
the 2k — 1 edges in some order. Letting T~ k be the finite set of shapes, R k can thus 
be described by a shape t S T 2 * fc and the edge lengths x±, . . . , X2k-i > 0, and for the 
Brownian CRT, Rk has density [HI Lemma 21] 

2fe-l 

f(i;x 1 ,...,x 2 k-i) = se~ s/2 , s = ^ x { . (5.2) 

i=i 

Aldous [SI Corollary 22] also gives a construction of the Brownian CRT in terms 
of a (normalized) Brownian excursion B ex . Let £ = 2B ex . Then Aldous shows that 
there exists a function £ mapping [0, 1] onto the Brownian CRT, with the Lebesgue 
measure mapped to v and, cf. (|5,ll) . 

d{C(s), C(*)) = CM + C(t) " 2m(s, t; () = a 2 (s, t; (). (5.3) 

Indeed, the Brownian CRT can be defined as the quotient space of [0, 1] with the 
semi-metric a 2 (s,t; £), identifying points of distance 0, see |21| . 

The function £ is not injective, but if C(s) = C(^) 5 and thus a 2 (s,t;() = 0, then 
C( s ) = C(^) = m ( s ,t',C)i which implies W(s) = W(t). Hence we can define a contin- 
uous random function W on the Brownian CRT by W(((s)) = W(s); conditioned 
on the CRT, the W(x) are jointly Gaussian with mean and, by 1)5. 1J) and 1)5.3)1 . 
Var (W(x) — W(y)) = d(x, y). Thus W is the random mapping of the Brownian CRT 
into R considered by Aldous [4; Aldous defines ISE as the measure on R that v is 
mapped to by W. This is clearly the same as the measure that W maps Lebesgue 
measure on [0, 1] to, i.e. the occupation measure of W as claimed above. 

6. Existence of the density: Proof of Theorem 12.11 

Although Theorem 12.11 follows easily from Theorem 13.11 and its proof, we find it 
interesting to give a different, self-contained proof. We use the standard Fourier 
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method, see e.g. JH] and the references there, together with Aldous's theory of the 
Brownian CRT [21 Oj . We define the Fourier transform /i of a finite measure [i by 
: = / e itx dfi(x). 

Lemma 6.1. If < a < 3/2, then 



E 



/oo 
(itri/i ISE (t)i) 2 dt < oo. 
-oo 



Proof. Since // ISE is the occupation measure of W, the head of the Brownian snake, 
its Fourier transform can be expressed as 

/oo rl 
e itx dfi 1SE (x) = / e itw ^ds. 
-oo io 

Consequently, |ju ISE (i)| 2 = J* /„* e i*(^( s )-^(«)) rf s d«. Conditioned on C, W(s)-W F («) 
is by (j5.1|) a Gaussian random variable with mean and variance o~ 2 (s, u; (). Hence, 

E(|/I ISE (t)| 2 | C) = C (\(e*mu)-W{s)) | ^ dsdu 



JO 

^^( S ,«;C)/2 ds ^ 



'0 Jo 

and thus, letting U\ and C/2 be independent uniform random variables on [0, 1], 

E\%^{t)\ 2 =E J 1 [\- t2 ^^y 2 dsdu = Ee- t2 ^ u ^y 2 . (6.1) 
io Jo 

Let C be as in Subsection 15.21 Then Xi := C(Ui), i = 1,2, are two independent 
random nodes in the Brownian CRT, and 1)5.3)) shows that 1)6.1(1 can be written 

E|/2 ISE (t)| 2 =E e -' 2 ^iA)/2. ( 6 .2) 

For a > we thus have, letting y = t 2 , 



E 



/"oo roc 

/ (t a \$ m (t)\)*dt = E t 2a e- t2d ^' x ^ 2 dt 
Jo Jo 

POO 

= E \y a - x l 2 e-^ X ^l 2 dy 



2t 

a-1/2. 



= Ei(d(X 1 ,X 2 )/2)- Q - i/ 'r(a + l/2) 
= C7i(a)Ed(X 1 ,X 2 )- a " 1 /2. ( 6 .3) 

From ()5.2() (with k = 2) follows the symmetry d(Xi,X<2) = d(Xi,o). Moreover, 
by the same formula 1)5. 2|) with A; = 1, d(Xi,o) has a Rayleigh distribution with 
density xe~ x I 2 . Hence, 

Ed(X 1 ,X 2 y a - 1 / 2 = Ed(X 1 ,o)- a -V 2 = / x- a -V 2 xe- x2/2 dx < 00, 

JO 

when a < 3/2, and the result follows from ()6.3)) and the symmetry of | /2 ISB | . □ 
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By Lemma l6~Tl if < a < 3/2, then (|y| a |/2 ISE (y)|) 2 dy < oo a.s. Taking first 
a = 0, we see that jj, iSE G L 2 (IR); by Plancherel's theorem |28| Theorem 7.9] this 
shows that /% E is absolutely continuous with a density / ISE G I? . Note that the 
Fourier transform / ISE coincides with j£t ISB . 

For a > 0, we define the (generalized) Sobolev space L 2,a by 

/OO 
((l + \t\r\f(t)\) 2 dt<oo}, (6.4) 
-OO 

where / is the Fourier transform of /. Lemma 16. II thus shows that a.s. / ISE G L 2 ' Q 
for every a < 3/2. (There is no problem with null sets, since it suffices to consider 
rational a, say.) 

Further, for < a < 1, we define the Holder space H a as the space of bounded 
continuous functions / on R such that \ f(x) — f(y)\ < C\x — y\ a for some C and all 
x and y. 

To show that / ISE is (i.e. can be chosen) continuous with the regularity properties 
in Theorem 12. 11 we use some general embedding properties of these spaces. 

Lemma 6.2. (i) If < a < 1/2 and 1/2 > 1/p > 1/2 - a, then L 2 ' Q C L p . 

(ii) // 1/2 < a < 3/2, then L 2 <° C H a _ 1/2 . 

(iii) If a > 1 and / G L 2 ' Q , £/ien / has a derivative f in distribution sense and 
a.e., with f G L 2 ' Q_1 . 

This lemma is well-known: (i) and (ii) are special cases of the Sobolev (or Besov) 
embedding theorem, see e.g. Theorem 6.5.1] or |291 Chapter V]; indeed, we may 
also take 1/p = 1/2 — a in (i). However, since the proof of the general embedding 
theorem is quite technical, we give a simple proof of this special case. 

Proof, (i): We may assume p > 2 since the case p = 2 follows by Plancherel's 
theorem. Define p' G (1,2] by 1/p' = 1 — 1/p. By Holder's inequality, 

/oo . poo \ p' /2 

\f\ p '<( {{i + \t\r'\f{t)ff /p dt) 
-OO i/ — OO ' 

x(r((i + iti)-^) 2/(2 ^' ) ^) lv/2 

J — OO ' 

= ii/iiC4/ji + i*i)- 2V/(2 " p ' ) ^ 1 ~ p72 <-' 

since it is easy to check that 2ap' > 2 — p' when 1/p > 1/2 — a. Consequently, 

/ G L p ' , which by the Hausdorff- Young inequality yields / G L p . By the inversion 
theorem for the Fourier transform (defined for tempered distributions, say), this 
yields / G IP . 

(ii): First, by Holder's (Cauchy-Schwarz's) inequality, 

/oo . /'oo . ]/2 / poo . ]/2 

i/i<(/ ((i + \t\r\f(t)\) 2 dt) (/ (i + |t|)- 2 ^) <oo, 
-oo J — OO J — OO 

since 2a > 1. Hence / has an absolutely integrable Fourier transform, which shows 
that / is a continuous bounded function given by the inversion formula f(x) = 
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(27r) 1 f e lxt f(t) dt. Hence, for any x and h > 0, using Holder's inequality again, 

1 f°° 

\f(x + h)-f(x)\ = — y_ (e-^+W-e-^fftdt 



oo 



< ( / i e ^-i| 2 iti- 2 «^y /2 f r \t\ 2a \f(t)\ 2 dt' 1 ' 



: (/I 2 "" 1 / \e lu - l? \u\- 2a du 



■oo 
1/2 

Wfha 

<C 2 (a)\\f\\ 2 , a h a -^ 2 . 

(iii): We have f'(t) = —it/, with /' taken as a distribution. Since / G L 2 ' a , this 
shows that /' G L 2 and thus /' € 1? by Plancherel's theorem. Consequently, by ele- 
mentary distribution theory, the derivative exists a.e., and equals the distributional 
derivative /'. Further, from the definition (|6.4[) . /' G L 2,a_1 . □ 

Since, as remarked above, Lemma 16 . 1 1 shows that a.s. / ISE G L 2,a for every a < 
3/2, Theorem 12. lT ii^ follows by Lemma I6.2f iih while Theorem 12 . lT iii) follows by 
Lemma E2f iii) and (i) (applied to // SE ). 

Finally, Theorem 12. 1H ) follows because ^ ISE has compact support, viz. the image 
of the compact set [0, 1] by the continuous function W. □ 

7. Local limit law for the density: Proof of Theorem 13.11 

The proof is based on the known convergence (jl.lj) of random measures. To obtain 
the stronger result in Theorem l3.1l on convergence of densities, we use a compactness 
argument as follows. We begin with a measure-theoretic lemma. Recall that a Polish 
space is a space with a topology that can be defined by a complete separable metric. 
For generalities on convergence of random elements of metric spaces (equipped with 
their Borel c-fields), see e.g. Billingsley |^j or Kallenberg ^B]- in particular, recall 
that a sequence (W n ) of random variables in a metric space S is tight if for every 
e > 0, there exists a compact subset K C S such that P(H / n G K) > 1 — e for every 
n; in a Polish space, this is equivalent to relative compactness (of the corresponding 
distributions) by Prohorov's theorem [Si Theorems 6.1 and 6.2], |18| Theorem 16.3]. 
Recall further that "convergence in distribution" really means convergence of the 
corresponding distributions, but it is often convenient to talk about random variables 
instead of their distributions. 

Lemma 7.1. Let S± and S2 be two Polish spaces, and let <j> : S\ — > £2 be an injective 
continuous map. If {W n ) is a tight sequence of random elements of S± such that 

<j>(W n ) — —> Z in S2 for some random Z G S2, then W n — —> W in S\ for some W 

with cf)(W) = Z. 

Proof. By Prohorov's theorem, each subsequence of (W n ) has a subsequence that 
converges in distribution to some limit. Let W' and W" be limits in distribution of 

two such subsequences W n > and W n >/. Since <fi is continuous, (j)(W n >.) <f>(W) and 
4>(W n ») <t>(W"). Hence, <f>(W) = Z = <f>(W"). 
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Let A be a (Borel) measurable subset of S\. By the Souslin-Lusin theorem [HI 
Theorem III. 21, see also III. 16-17], 4>{A) C £2 is measurable. Thus, using the 
injectivity of 0, 

F(W' G A) = F((j)(W') G <p(A)) = F(<P(W") G <p(A)) = F(W" G A). 

Consequently, W = W" . 

In other words, there is a unique distribution of the subsequence limits. Thus, if W 
is one such limit, then every subsequence of (W n ) has a subsequence that converges 

in distribution to W; this is equivalent to W n — > W . □ 

Let Y„ denote the random probability measure on the left-hand side of 1)1. ljl . let 
be the probability measure with the triangular density function h~ l {l — \x\/h)+, 
and let Y n be the convolution Y n * i/ -1/4. Note that y n has the density g n (x) := 

?i- 1 7~ 1 n 1 / 4 X n ( 7 - 1 n 1 /4 x ) G C (R). Since Y n MlSB by flDJ, and V - 1/4 <5 , 

it follows easily that Y n fi 1SE too. 

Let <Si := {/ G Co(K) : / > 0}, with the uniform topology inherited from Cq(R), 
let ^2 be the space of locally finite measures on R with the vague topology, see e.g. 
Kallenberg |18l Appendix A2], and let 4> map a function / to the corresponding 
measure f dx, i.e., <f>{f) is the measure with density /. Then S\ is a closed subset of 
the separable Banach space Co(R), and is thus Polish, and so is £2 by |18l Theorem 
A2.3]. Further, <fi is continuous and injective. 

Take W n := g n in Lemma 17.11 We have just shown that 4>{g n ) = Y n /x ISB in 
the space of probability measures on R and thus also in the larger space ^2 . If we can 
show that the sequence g n is tight in Si, or, equivalently, in Co(R), then Lemma l7.ll 

shows that g n — —> g for some random function g G Co(R), which further equals (in 
distribution) the density / ISE of Hise] hence the conclusion of Theorem 13. II follows. It 
thus remains only to prove the following lemma. 

Lemma 7.2. The sequence g n (%) ■= n~ 1 7~ 1 n 1 / 4 A > n (7~ 1 n 1 / 4 x), n = 1,2,..., is 
tight in C (R). 

The central estimate in the proof of Lemma 17.21 and thus of Theorem 13. 1| is the 
following, which will be proved in Section |H1 For a sequence x(j), we define its 
Fourier transform by x(u) := x(j)e 1JU ; this equals the Fourier transform of the 
measure ^ ■ x(j)5j on R. 

Lemma 7.3. There exists a constant C\ such that for all n > 1 and u G [— it, it], 

V\n- X Xn{u)\ 2 <-^- A . (7.1) 
1 + nr 

We can now prove Lemma 17.21 as follows. We have Y n (y) = n~ 1 X n (7re -1 / 4 y) . 

Consequently, Y n is a periodic function with period 2-7T7~ 1 n 1 / 4 , and Lemma 17.31 
translates to 

E |^(y)| 2 < t-^Vt, \V\ < 7- V /4 vr. (7.2) 

1 + 7 V 
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Further, 

g n {y) = %{y) = ^(y)v-v4(y)- (7.3) 

Lemma 7.4. Suppose that < a < 3. Then there exists a constant C(a) such that 
if h > and f is a function with period 27r/h, then 

\y\ a \f(y)\ 2 \My)\ 2 dy < C(a) [** \y\ a \f(y)\ 2 dy. 

J-w/h 

Proof. By the change of variables y i— » h~ l y, we may assume h = 1. Then Pi(y) = 
(sin(y/2)/(y/2)) 2 . Hence, for / and |y| < tt, 

r ^ot. m sin2 (^/ 2 ) ^ y 2 



and 



(2fc+l)7T /*7T 4 

b| a |/(y)| 2 pi(y)| 2 <fo < (3n\k\) a / l/d/)! 2 ^ dy 

<Ci(o)|fcr 4 T M a |/(l/)| 2 dy. 

For the case k = 0, we use instead the estimate |^i(y)| < 1. The result follows by 
summing over all /c. □ 

Let h : = 771 -1 / 4 . Then, by Lemma 17.41 (|7.3|) and (|7.2j) . for any fixed o with 
< a < 3, 

/"OO P7r/h 

E \y\ a \g n {y)\ 2 dy<C{a)E \y\ a \Y n (y)\ 2 dy 



n/h 
w/h 



= C(a) r \y\ a E\Y n (y)\ 2 dy 
J-n/h 



We have proved the following, taking a = 2a. 

Lemma 7.5. If < a < 3/2, then E ||s„||| a < C(a), for some C(a) not depending 
on n. □ 

Next, fix p G (0,1), and let a = f3 + 1/2 < 3/2. For A,Af > 0, let if M ,A be 
the set of all functions / in Cq(M) such that f(x) = for \x\ > M and ||/||2,a — 
A. By Lemma l6.2f iiL the functions in Km, A are all Holder (/3)-continuous with 
uniformly bounded norm; they thus form an equicontinuous family. We may regard 
Km,a as a subset of C[—M,M], the space of continuous functions on the compact 
interval [—M,M], and it follows by the Arzela-Ascoli theorem |28| A5] that Km, A 
is a relatively compact subset of C[—M,M], and thus of Co(M) too. (Note that the 
functions in Km, A all vanish at ±M.) 



16 MIREILLE BOUSQUET-MELOU AND SVANTE JANSON 

Let e > 0. It follows from Lemma 17.51 that there exists A such that IP(||gn||2,a > 
A) < e/2 for every n. Moreover, Marckert |231 Theorem 5] also showed that 

n- 1/4 sup{|i| : X n (j) + 0} = n" 1 / 4 sup{\£(v)\ : v E T n } -±> W (7.4) 
for some random variable W. It follows from (|7.4|) that there exists M such that 

¥(g n (x) 7^ for some x with \x\ > M) 

= P(X n (j) / for some j with \j\ > 7~ 1 n 1/4 M - l) < e/2. 

Consequently, P(<? n 6 Km,a) > 1 — £ for every n, which shows that the sequence (g n ) 
is tight. This completes the proof of Lemma 17.21 and thus of Theorem 13.11 except 
for the proof of Lemma 17.31 

Remark 7.6. A more concrete alternative to the compactness argument (Lemma l7.1|) 

used above is to define regularizations of functions / on R by f^ h \x) := /i" 1 f 

for h > 0. Note that fj^ '(x) = /i _1 /i ISE [x, x + h]. For each fixed h > 0, the version of 
(|3.1j) with both left- and right-hand side regularized holds since the corresponding 
distribution functions converge. We may then let h — > 0, using the Holder estimate 
obtained by Lemmas 17.51 and 16.21 together with Theorem 4.2]. 

8. Proof of Lemma 1731 

It remains to prove Lemma 17.31 We consider first the case of binary trees. We 
introduce the sequence of generating functions 

TeT i=l v&T 

where T is the family of all (possibly empty) binary trees and \T\ is the number of 
nodes in T. Thus is a power series in t, with coefficients in Z[xi, . . . , x^, l/x±, . . . , 1/x/.] 
the ring of Laurent polynomials in the Xi with integer coefficients. For k = 0, the 
product in the definition of Fq reduces to 1, so that Fq is simply the generating 
function of binary trees. In what follows, we often denote x = (x\, . . . ,Xf-) and 
Fk{x) = Fk(t, xi, . . . , Xk). Moreover, for any subset I of [k] := {1,2, ... ,k}, we 
denote xj = (xi t , Xj 2 , . . . ,X;) if / = . . . , i p } with i\ < • • • < i p . 

Proposition 8.1. The series F^ can be determined by induction on k > using 

F k (x) = l [k=0] +tJ2 (n**J \ U x n F\i\(xi)F\J\(xj) (8-2) 
(j,j) Kiel J \jeJ J 

where the sum runs over all ordered pairs (/, J) of subsets of [k] such that If] J = 0, 
and Xi = l/x{. In particular, 

1 - jT^M 



2t 

and each F k (x) admits a rational expression in Fq and the X\ 



1.3) 
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Proof. The equation satisfied by Fq reads Fq = 1+tF^ and is of course very classical: 
it is obtained by splitting a binary tree into its left and right subtrees. Note that the 
empty binary tree does not contribute to Fk when k > 0. Then, every non-empty 
binary tree is formed of a root with a left subtree T\ and a right subtree T 2 . Hence, 
for k > 1, 

Fk (x) = £ n ( 1 + e *f M + e 

Ti,T 2 i=l \ ueTi KgT 2 



E^ l+Pil EII E^" 1 

Ti,T 2 (J, J) i€/ VdGTi 




where the sets / and J are as in the statement of the proposition. The result follows 
upon exchanging the two sums. □ 

Actually, for the proof of Lemma l7.Hl we need only a special case of F2. The above 
proposition gives a simple explicit expression of i*2(i, x, y) in terms of Fq (and x and 
y), or, equivalently, in terms of the generating function 5 = Fq — 1 of non-empty 
binary trees: 

fl-aw- 1 -"-/ 3 * . ,8.4) 



F 2 (t, e '",e-'")= ^ „T„ , „ op ...A2 - < 8 ' 5 ) 



Corollary 8.2. For any real it, 

5(l + 5)(l + 25-5 
(1 - B)(l + B- 25 cos it 

Proof. The cases fc = 1 and = 2 of the previous proposition give 

Fi(z) = iF 2 + t(x + x)F Fi(x) 

and 

F 2 (x,y) = iF 2 + t(x + x)F F!(x) + % + y)F F 1 (y) 

+ + xy)F 1 (x)F 1 (y) + i(xy + xy)F F 2 {x,y). 

Using F = 1 + B and i = 5/(1 + B) 2 , this gives 

5(1 + 5) 



F x (x) 

5 2 (x,y) 



l + B(l-x-x)' 
B(l + B)(l + 2B + B 2 (l -xy- xy)) 



(1 + B (l - x - x))(l + 5(1 - y - y))(l + 5(1 - xy - xy)) ' 
Specializing to x = l/y = e lu provides the result. □ 
By definition, 

X n {u) = Y J X{T,T n )e^ = £ e M M«. (8.6) 

j veT n 
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Hence, if T n := {T £ T : \T\ = n} is the family of binary trees of size n, 

TeT n v£T n 

and 



F 2 (t,e iu ,e- iu ) = J2t n \T n \K\X n (u)\ 2 . (8.7) 



n=l 



Since \T n \ = [t n ]B(t) = ^ T ( 2 n ") ~ 7r _1 / 2 n- 3 / 2 4 n , fCB is equivalent to 

n i/2 

n^tM 1 "^ 1 ") < C 2 4™-— -j, | U | < tt. (8.8) 
We will prove this using complex analysis. We begin by studying B. 

Lemma 8.3. B = B(t) is a bounded analytic function of t in the domain D : = 
C \ [1/4, +oo). Moreover, for t 6 &D = [1/4, +oo), B has continuous boundary 
values B + {t) and £?_(i) from the upper and lower side. B(t) (extended by B + or 
B-) is real if and only if t € (— oo, 1/4]; on this interval B(t) is strictly increasing 
from -1 to 5(1/4) = 1. 

Proof. The first assertions are immediate from (|8.4|) . Next, if B(t) is real, then so 
is t = B/(l + B) 2 . It follows further from (jHHJ) that B(t) is real for t < 1/4, but 
B±(t) is not real for t > 1/4. The formula t = B/(l + B) 2 shows further that 
B = — 1 is impossible, and that B = 1 if and only \i t = 1/4. Since B(t) — ► —1 as 
t — ► — oo, it follows by continuity that —1 < B(t) < 1 for t < 1/4. For such t we 
have d-B/di = (dt/dB)' 1 = (1 + -B) 3 /(l - 5) > 0, which completes the proof. □ 

Let us for simplicity write F u (t) := Fi(t, e lu , e~~ xu ). 

We first observe that, for any real u, F u (t) is an analytic function of t in the 
domain D' := D \ (— oo, — 3/4]. Indeed, by Corollary 18.21 and Lemma 18.31 F u (t) is 
meromorphic in T> with poles when 1 — _B = 0orl + J3 — 2B cos u = 0. In the first 
case, B = 1 and thus t = B(l + B)~ 2 = 1/4, which is outside V. In the second case, 
B = 1/(2 cos u — 1). Since 2cosu — 1 E [—3, 1], this means that B is real and either 
B > 1 or B < —1/3. By Lemma 18.31 B > 1 is impossible in P, while -B < —1/3 
implies t = B(l + B)- 2 < -3/4. 

At this stage, we can apply, for any fixed value of u, the standard results of 
singularity analysis J2]. For u = 0, we find 



4 n n i/2 

[t n ]F u (t) = —=- (1 + 0(l/n)) 

V7T 



while for 



These results are certainly compatible with the desired bound (|8.8j) . but, as we need a 
uniform bound, valid for all u, we have to resort to the basic principles of singularity 
analysis. 
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By the Cauchy integral formula, 

for any contour T in T>' that loops once around 0. We assume n > 4 and choose a 
contour T = rW that depends on n: T^ n ' := T± UT + UT2 UT_, where Ti is the circle 
\t — 1/4 1 = l/n in negative direction, T + and T_ are both the interval [1/4+1/n, 1/2], 
taken in opposite directions and using the boundary values B + and S_, respectively, 
and T2 is the circle |t| = 1/2 in positive direction. (For convenience, we have pushed 
the contour to include part of the boundary of T>'; the reader who prefers staying 
strictly inside T>' may replace T± by line segments close to the real axis.) 
Next we estimate \F u (t)\ on r( n ). 

Lemma 8.4. For all t G r( n ) and u G [— ir, it], 

n 3/2 

\F u (t)\<C 3 — — (8.10) 
1 + nir 

Proof. We claim that for t G rW 

|1 - B{t)\ > cin~ 1/2 , (8.11) 

1 1 + B(t) — 2 cos uB(t) I > C2 max(n" 1/2 ,l-cosn). (8.12) 

The result then follows from 1)8. 5 Jl and 1 — cosn > c^u 2 . 

In fact, (|8.11|) is the special case u = of (|8.12() . so it suffices to prove the latter. 
Since, by compactness, \B(t)\ > C4 > for t G T C {t : ^ < \t\ < ^}, it is enough to 
prove 



\Bity 1 + 1 - 2 cosn) > c 5 max (n" 1/2 ,l- cosn). (8.13) 



Indeed, 



5 t " 1 + 1 - 2 cosn = - + 1 - 2cosn 



l ~ At+ 2 f^ + 2(1 -cos,). (8.14) 



For £ G Ti, this is 2^1 - 4t + 0(l/n) + 2(1 - cosn). Since Re y/1 - At > and 
1 — cos n > 0, then 

15(f)- 1 + 1 - 2cosn| > |2 v / l-4t + 2(l-cosn)| -0(l/n) 



> max{2| v / l - 4t|, 2(1 - cosn)} - 0(l/n) 
= max{4n~ 1/2 ,2(l - cosn)} - 0(l/n), 

which yields gggj) . 

On T-i-, \/l — 4t is imaginary, and 

llm^C*)- 1 + 1 - 2 cos n) I = \^±E3. > ^li^l > 2r l - 1 / 2 . 



.15) 
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Moreover, if further 1 - cos u < 2^4t - 1, lj%T5j) also yields |lm(.B(i) _1 + 1 - 
2 cos uj | > y/4t — 1 > ^(l ~~ cosu). If, on the contrary, 1 — cosu > 2\/4t — 1, 
then, because < 4t — 1 < 1 and thus 4t — 1 < \/4i — 1, 

|Re(S(t) _1 + 1 - 2 cos u) I = 2(1 - cosu) - ^ 

> 2(1 - cosu) - 2{4t - 1) > 1 - cosu. 

In both cases, l-B(i)" 1 + 1 — 2 cos u\ > ^(1 — cosu), which together with (|8.15|) 
completes the verification of ()8.13j) for t G T±. 

Finally, for t G I^, we use compactness. We observed above that 1 + B(t)(l — 
2cosu) = is possible only for t G ( — 00,— 3/4] U {1/4}, and in particular not for 
t G hence fnf ier2 >ug [o,27r] \B{t)~ 1 + 1 — 2cosu| = cq > 0, which implies (j8.13j) for 
t G IV □ 

By and (l8~TUl) . 

WW < / I^WI ir"" 1 1^1 < £37^—4 / Itr"" 1 |««|. (8.16) 

For t G Ti, It]-™- 1 = 0(4 n ), and thus J Fi \dt\ = 0{n~ 1 4 n ). Secondly, 

Jr± I*! - " -1 < !t/lt~ n ~ X dt < n~ l 4 n . Finally, /p 2 |t|- n - 4 \dt\ = 0{2 n ). Sum- 
ming these estimates we find / r \t\~ n ~ 1 \dt\ = 0(n _1 4 n ), which together with (IH7TH1) 
completes the proof of (|8.8j) and thus Lemma 17.31 in the case of binary trees. 

For complete binary trees, we use the well-known equivalence between binary and 
complete binary trees, where a binary tree T of order n is identified with the internal 
nodes in a complete binary tree T c of order 2n + 1. With this identification, one has 

(X(j -l;T)+X(j + l;T) ifj/0, 
[1 + X(-1;T) + X{1;T) otherwise. 

Hence, temporarily using X% instead of X n for the complete binary trees, it follows 
from (|8.6|) that X% n , 1 (u) = 1 + 2 cos uX n (u). Hence the estimate in Lemma f7.3l holds 
for complete binary trees too (possibly with a different constant). 

9. Proof of Corollaries I3.2hl3.4l 

Proof of Corollary 13.21 This is immediate from (|3.2(l and the fact that f n — > f in 
Co(M) implies f n {jn/n l ^ A ) — ► f(x), see e.g. Theorem 5.5]. □ 

Proof of Corollary 13.31 By symmetry, fi SE (x) = f 1SE (—x), so we may suppose x > 0. 
Then, as shown by Bousquet-Melou §6.2.2, Conjecture 15 and Theorem 14], for 

naturally embedded random binary trees, X n ([xn 1 / 4 J )/n 3 / 4 — > 2~ 1 / 2 y(2~ 1 / 2 x) 
for a family of random variables Y(u), u > 0, with moment generating functions 

Ee° y M = L(u,a). Combining this with Corollary EHl we find 2- 1 / 4 / ISE (2- 1 / 4 x) = 
2" 1 / 2 y(2- 1 /2 x ) ) x > 0, and thus /i SE (x) = 2- 1 / 4 y(2" 1 / 4 |x|). (The normalization of 
/ise in Hi is different.) □ 



X(j;T c 
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Proof of Corollary EM By the proof of Corollary EH / ISE (0) = 2~ 1 / 4 Y(0), where, 

by Proposition 12 and Theorem 14] Y(0) = 2 1 / 2 3~ 1 T~ 1 / 2 . 

The (negative) moments of T are given by the standard formula E T~ s = T(3s/2 + 
l)/r(s + l), s > -2/3. □ 



10. Other tree models 

Consider a randomly labelled conditioned Galton- Watson tree as in Coniecture l3.5l 
We know that the global limit result holds, and the proof in Section holds ver- 
batim in this case too and shows that to prove Conjecture [Hill it is sufficient to verify 
that the estimate of Lemma 17.31 holds. We have not been able to do so in general, 
but we can show the required estimate in the two special cases in Theorem 13.61 

We consider thus in this section the two families of labelled plane trees that were 
studied in [7j. In the first family T 1 , the root is labelled 0, and the labels of two adja- 
cent nodes differ by ±1. In the second family T 2 , the latter condition is generalized 
by allowing the increments along edges to be 0, ±1. 

Again, we introduce a sequence of generating functions: 

k 

F k (t, Xl ,..., x k ) = F(x) := £ t^ H(E x f V) ) > ( 1(U ) 

TgT i=l veT 

where T is either T 1 or T 2 and \T\ is the number of edges in T. The following 
proposition is the counterpart, for each of the two new families, of Proposition 18. II 

Proposition 10.1. For plane trees with increments ±1, the series Ff. can be deter- 
mined by induction on k > using 

F k (x)=i+tY, (n> + n x *) f \i\( x i) f \j\( x j) 

IC[k] i€l / 

where J = [k] \ I and = 1/xj. For trees with increments 0, ±1, the above equation 
becomes 

F k {x) = l + tJ2 1 1 + Y[xi + Y[xi) F w { Xl )F m (xj), 
ic[k] V iei iei ) 

with the same notation as above. In both cases, each F k {x) admits a rational expres- 
sion in Fq and the Xi. 

Proof. The proof is very similar to that of Proposition l8.ll We now use the standard 
recursive description of plane trees based on the deletion of the leftmost subtree T\ of 
a tree T (not reduced to a single node). This leaves another plane tree T%. Also, one 
has to take into account the fact that the label of the root of T\ may now take two 
(or three) different values (depending on the family of trees under consideration). 
Finally, the tree reduced to a single node contributes 1 in each F k . When T = T 1 
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and k > 0, this gives 

F k{x) = i + y: * 1+|Tll+|T2 ' n ( E + E x 

Ti,T 2 t=l 

+ ^^i™niE^ +1 +E-? 

Ti,T 2 i=l »eT 2 

and the result follows after expanding the products, and then exchanging the sums. 

□ 

We easily find explicit formulas for F\ and F2 from Proposition 110, ll cf. the proof 
of Corollary 18,21 We leave the details to the reader and state only the result that 
we need, in terms of the series T = T(t) that counts labelled trees not reduced to a 
single node. Depending on which tree family is studied, one has 

T = TW := B(2t) = j^jWjZg for r i 

v ; At 

T = TV :=B m = 1 - 6t -^ Wt forT 2 , 

6t 

where the series B(t) is defined by 1)8.4(1 . 

Corollary 10.2. For plane trees with increments ±1, 

F2(t ' 6 ' 6 } " (l-T)(l-Tcos^ • (10 ' 2) 
For plane trees with increments 0, ±1, 

F2(t ' 6 ' 6 } " (l-T)(3-T(l + 2cosu)) 2 " (10 ' 3) 

We may now complete the proof of Theorem 13.61 by the argument in Section |HJ 
we give a sketch only and leave again the details to the reader. First, the functions 
F 2 (t/2,e iu ,e- iu ) (for T 1 ) and F 2 (t/3, e iu , e~ iu ) (for T 2 ) are analytic functions of 
t E D for every real u. Next, in analogy with Lemma 18.41 with the same contour 
r( n ) as there, for t G T and \u\ < vr, 

\F 2 (t/2,<r,e- m )\ < C A -—- A for T , 

1 1 1 + nu 4 

„3/2 

|F 2 (V3, e m , e - m )| < C 5 — - — for T 2 . 

1 1 1 + nu 4 

Indeed, the proof is almost exactly the same; we replace the left-hand side of 1)8.12(1 
by 1 1 — cos uB(t) I and 3 — B(t) — 2 cos uB(t) | and similarly the left-hand side 1)8.13(1 
by \B(t)~ l — costi| and |3-B(i) -1 — 1 — 2 cos ii , note the corresponding changes in 
1)8.15)) and argue as before. 



DENSITY OF THE ISE AND LIMIT LAWS FOR EMBEDDED TREES 



23 



The Cauchy integral formula (|8.9j) then leads to, cf. (|8.8j) . for |u| < ir, 

n 1 / 2 

[t n ]F 2 (t,e lu ,e-™) < C 6 8 n -—^, for T 1 , 

1 + nu^ 

[t n ]F 2 (t,e iu ,e- iu ) < C 7 12 n -^^, for T 2 . 

1 + nnr 

By (jST) and \T^\ = 2 n [t n ]B(t) ~ 7 r- 1 / 2 n- 3 / 2 8 n , |T n 2 | = 3 n [t n ] J B(t) ~ Tr -1 / 2 ^ 3 / 2 ^", 
this yields i|7.1j) for these two families. (Note that we have let \T\ be the number 
of edges for T 1 and T 2 ; thus we now should replace X n by X n+ \ in (|8,7j) . but this 
makes no difference for (|7.1jl .) 

This completes the proof of Theorem 13.61 



11. Moments of the density of ISE 

We know by Corollary 13 .31 that /ise(A) has a moment generating function (defined 
in an interval containing 0), and thus finite moments of all orders. We next present a 
formula for these moments, and more generally for mixed moments involving several 
values of A. We use a general method for occupation densities of Gaussian processes. 
To state the formula, we introduce more notation. 

Given Q (which as always is 2B CK ), and k points si,...,s k £ [0,1], the random 
vector (W(s\), . . . , W(sh)) has a Gaussian distribution with mean and covariance 
matrix 

S C;Slv .., Sfc := {m{s u sf C))*j =r (H- 1 ) 

We let (f( ]Sl ,...,s k denote the density function of this distribution. (We may ignore 
the cases when the distribution is degenerate; a.s. this happens only when Sj = Sj 
for some i and j.) 

Using the construction in Subsection 15.21 of the Brownian CRT from £, we can 
transfer these notations to the CRT. Given £ and k points x\,... ,x k in the corre- 
sponding CRT, the random vector {W{x\), . . . , W(xk)) has a Gaussian distribution 
with mean and covariance matrix 

Sc;x lv ..,x fc := {m{xi,Xj-X)) k itj=v (H-2) 

where m(x, y; Q is the length of the common part of the paths from the root to x and 
y in the CRT. We let <p^- xi ,...,x k denote the density function of this distribution, and 
note that if Xj = C( s «)j i = 1, • • • , fc, then m(xi,Xj] C) = m(si, Sj] () and ip£- lXl ,...,x k = 
¥>C;«i,.«,«fc- 

We further let X\, . . . ,X k denote k independent random nodes in the Brownian 
CRT. 

Theorem 11.1. For any real numbers X\, . . . , A^, 

E(/ ISE (Ai) • • • /i SE (Ajfc)) = E / ••• / (p ( . slj _ iSk (Xi, . . . , X k ) dsi ■ ■ ■ ds k 

Jo Jo 

= E<p c . Xl ^ Xk (X 1 ,...,X k ). (11.3) 
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Proof. The equality of the last two expressions follows by the construction of the 
Brownian CRT and the definitions above. 
We define, for A G R and h > 0, 



(11.4) 



Z h (X) = jf 1 [w(sMXA+h]] ds = fc-WlA, A + fc] 

= ^i" 1 fws(y) dy. 

Since / ISE is continuous by Theorem 12.11 Z^{X) — > /i SE (A) a.s. as h — > 0. From this 
definition follows 

E(Z h (Ai)---Z fc (A fc ) | C) 

••• / h~ k F(W(s t ) G [Ai,Ai + /i], i= l,...,fc | C)<fei— d»* 
•/ o 

1 /-l rXi+h rX k +h 

•••/ /i~ fe / •••/ tpt; S1 ,...,s k (yu---,yk)dyi--- dy k ds!--- ds k 

JO JAj JA fc 

Ai+h /■A fc +/i 



E /» 



Ai 



A* 



<Pc-,xi,...,x k (yi, ■ ■ ■ > 2/fc) dyi • • • dyfc 



and thus 



<PC;Xi,...,x k {yi,---,yk)dyi--- dy k . (11.5) 



E(Z fc (Ai) • ■ • Z ft (A fc )) 

rXi+h rXk+h 

= Eh~ k •■■ / 

JAi JA fc 

To obtain the conclusion, we now let /i — ► 0; however, we have to justify taking 
the limit inside the expectations on both sides. For the right-hand side, we use the 
fact that a Gaussian distribution in R fc with mean has a density function that has 
its maximum at 0; hence we can, by Lemma ll 1.21 b elow . use dominated convergence 
with ■ ■ ■ ,0) as dominating function. Since ip^-Xx,...,x k i s continuous, the 

right-hand side of (|11.5|) thus converges to the right-hand side of (|11.3j) . 

For the left-hand side we begin by applying Fatou's lemma, which now shows that 
the left-hand side of (|11.3|) is at most equal to the right-hand side. By Lemma 111. 21 
below, this yields a uniform bound, C k say, of the left-hand side for all Ai, . . . , A^. It 
follows from (tTOl) that E(Z h (A 1 ) • • • Z h (\ k )) < C k too, for every h > 0. If we here 
replace k by 2k, repeating every Aj twice, we see that the random variables Vh := 
Z/i(Ai) • • • Zh(X k ) satisfy E < C2 k - The variables Vh are thus uniformly integrable, 



follows EV h -> E(/ ISE (Ai 



/iSE(Afc)), 

□ 



and from V h ->• /i SB (Ai) • • • /i SE (Afc) as h 
see e.g. |14| Theorems 5.4.2 and 5.5.2]. 

Lemma 11.2. For every k > 1, M(p^.x lt ... x fc (0, • • • ,0) < oo. 

Proof. The subtree R k of the Brownian CRT spanned by X\ , X k and the root o 
has k — 1 internal nodes. Let R' k be the subtree spanned by o and the internal nodes 
of R k , and let l\, . . . , l k be the lengths of the k edges that attach X\,. . . , X k to R' k . 
The values of W along R k form a branching Brownian motion, i.e., W is a Brownian 
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motion along each edge of Rk and all increments are independent. In particular, 
conditioned on and the values of W on R' k , the values W(X±), . . . , W(Xk) at the 
leaves are independent Gaussian variables with some means b\ , . . . , and variances 
£1, ... ,£k- The conditional density function is thus at most W^(2itli)~ l l 2 , and thus, 



taking the expectation and using (15.2(1 

k 

E^i,-,x*(0,...,0)<EjJ(2 7 r£ i 



k 

1/2 



(2fc-3)!!(2vr)- fc /2 /" . . . /" \{^ 2 l^-^' 1 ^ 



2 



dl x ■ ■ ■ di 



2k-l 



< oo. 



□ 



Since the distribution of the covariance matrix S^ ; Xi,...,X/t is given by (|11.2|) and 
([5.2)1 . it is in principle possible to write the right-hand side of (|11.3|) as a multi- 
ple integral. However, the expression becomes rather complicated for higher mo- 
ments. In the simplest case Ai = • • • = A& = 0, (|11.3|) reduces to E/i SE (0) fc = 

(2-7r)~ fe / 2 E(det(S ( j.x 1 ,...,x fe ) 1 ^ 2 ), but even this seem difficult to compute in general. 
(These moments were found by another method in Corollary 13.41 ) 

In the case k = 1, Theorem 111.11 yields a simple formula for the average E/ ISE 
of the density, which equals the density of the average E/x ISE , i.e. the density of a 
random point chosen according to the random ISE. In the latter formulation, it was 
found by Aldous [I]. 

Corollary 11.3. For any real X, 

E/ ISE (A) = (27T)- 1 / 2 jf 2/ 1/2 exp(-^ - \) dy. 

Proof. <££ ; Xi(A) = {2-Ky)~ l / 2 e~ x2 /( 2y \ where y = d(X\,o), and y has the density 
function ye~ y2/2 by (JOJ. □ 

Alternatively, expanding the Laplace transform of Corollary 13.31 in a gives (see j7j 
Proposition 13]): 



E/ ISE (A) = — ^ ^ cos r 



m 



7r * — ' ml 

Both expressions yield E/ ISE (0) = 2~ 3 / 4 7r~ 1 / 2 r(3/4), as given by Corollary 13.41 

From Corollary 111.31 follows easily by integration another formula by Aldous [1]; 
we leave the proof to the reader. 

Corollary 11.4. For every real a > — 1, 

23a/4 



f°° f°° 2 3a/4 /a 1\ /a 

Ey \ x \ a dfi 1SE (x) = E J \ x \ a f 1SE (x)dx = _r(- + -)r(- + i 



□ 
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This extends (|4.5jl . (To see the equivalence when a = 2k, use the duplication 
formula for the Gamma function twice.) 

12. The grand-moments of the ISE: Proofs 

Proof of Theorem 14. H Let / be a continuous function on R. First, if / is bounded, 
then J f d[i is a continuous functional on the space of probability measures on 

R, and since [i n —> /i ISB in this space, see it follows that 

J fdfl n -±> y /d/W (12.1) 

We need to extend this to unbounded /. Thus, let f u , for u > 0, be the function 
that is equal to / on [—u, u], and is constant on (— oo, —u] and on [u, oo). Since f u is 

bounded, (|12.1j) applies to f u , i.e. J f u d[i n f f u d[i 1SE for every u > 0. Moreover, 
let V n ■= sup{\x\ : x £ supp/x n } = (2n) -1 / 4 sup{|£(u)| : u G T n }. By Marckert [23 
Theorem 5], V n V for some random variable V. Consequently, 

limsupP( / f u dn n + [ fdfin) < limsu P P(T/ n > u) < F{V > u), 

n— >oo \J J ' n— >oo 

which tends to as u — > oo. Finally, J f u dfi lSE — > f f dfi lSE as u — ► oo, since ^ ISE has 
compact support. Consequently, Theorem 4.2] shows that ([12. lj) holds for any 
continuous /. 

Taking f(x) = x k in 1)12.1(1 . we obtain the convergence m^^ n — ► of the mo- 
ments asserted in Theorem 14.11 Moreover, taking / to be a linear combination of 
such monomials, we see that joint convergence holds by the Cramer-Wold device 
Theorem 7.7]. 

In particular, for any partition A, m\ n — —> m\. We will show that the expectation 
n) converges as n — > oo. Applying this to the partition A' where each part in A 
is repeated twice, we see that also E(m| n ) = E(my )fl ) converges. The variables rn\ n 
are thus uniformly integrable, and the limit of their expectations E(m A n ) equals the 
expectation K(m\) of their limit, see e.g. |141 Theorems 5.4.2 and 5.5.2]. 

To complete the proof of Theorem 14.11 we thus have to show that the grand- 
moments E(m A n ) of /i n converge to the limits stated in the theorem. We introduce 
the non- normalized moments of fj, n : 

M*,n= J2 i ( V ) i = 2i/4nl+i/Arn i,n, 
v£T n 

as well as their factorial version, which is simpler to handle via generating functions: 
M hn = ^ *(v) Wv) -!)••• {i(v) -k + 1). 

v£T n 

We also use the notation M\ >n , analogous to (|4.2|1 . Then the remaining part of 
Theorem 14 . 1 1 easily follows from: 
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Proposition 12.1. As n — > oo, the non-normalized factorial moments of fi n satisfy 

E( ^" ) = r ( p + |A|/4-i/2) fe + °W> 

Proof. Let us first relate M\ <n to the generating functions of Proposition 18.11 It is 
simple to see that 

d x F p : = Xp (t, 1, . . . , 1) = £ t m f[ [ £ l {v)m - 1) • • • (/(„) 

OT X • • • OXp Te j- i=1 



■J&T 




n>0 

where C n = ( 2 ™) /(n + 1) is the number of binary trees with n nodes, known as the 
nth Catalan number. By Proposition 18. 11 the series d\F p is a rational function of t 
and \Jl — 4i. We want to study the singularities of these series. We will prove that, 
for p > 0, 

p (l-4t) e A ' v ' 

where P\(t) and Q\(t) are two Laurent polynomials in t, and 

e\=p 



1 


|A| 


1 fp + |A|/4- 


1/2 


if A is even 


2 


2 


2 | p +|A|/4- 


3/4 


if A is odd. 



(Note that P\ and Q\ may be singular at t = 0, although c^-Fp is analytic there.) 
From (|12.3|) . it follows that the only possible singularity of d\F p is at t = 1/4, and 
that, as t — ► 1/4, 

ca + o(l) 

A P (1 _ 4t)^+|A|/4-l/2 ' 

where c\ = P(l/4) when |A| is even and c\ = when |A| is odd. We will further 
show that the numbers c\ satisfy the recurrence relation (|4..3ft . The form (|12.3|) and 
the above singular behaviour do not hold when p = 0, and should be replaced in this 
case by the expression (|8.3|) of Fq and the singular behaviour 

F = 2 - 2y/l - At + 0(1 - it). 

Assume for the moment that we have proved (|12.3j) . Then the standard results of 
singularity analysis provide 

An p+|A|/4-3/2 

[t n ]d x F p = C n E(M x>n ) = r(p+ | A | /4 _ 1/2) (ca + o(l)) • 

Given that C n ~ 4 n n~ 3 / 2 /r(l/2), this gives the result stated in the proposition. 
Note that this asymptotic behaviour also holds for p = 0, with C0 = —2. 

Let us thus focus on (|12..3j) . Our proof works by induction onp+ |A|. 

• If p = 0, then A is the empty partition, and we have worked out above the value 
of Fq and its asymptotic behaviour when t — > 1/4. 
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If p > and Ai = 0, then 



d\F p = ^i"C re E(M A , n ) = £V<7 B nE(M A> ) = t^d y F p ^, 

n>0 n>0 



where A' = (A2, • • • , A p ). Then the form ()12.3|) follows by a simple calculation from 
the induction hypothesis, and the fact that e A = e\> + 1. (We do not give the details.) 
This calculation also provides the value of c A in terms of c\>. The case p = 1 and 
A = (0) has to be treated separately, since in that case A' = and the form 1)12.3(1 is 
not valid. 

• If p > and Ai > 0, then all the parts of A are positive. Let us differentiate 
()8.2j) Ai times with respect to x%, then A2 times with respect to X2, and so on, and 
then set Xi = 1 in the result. Since Aj > for all i, the terms for which I U J 7^ [p] 
do not contribute, and we are left with 

d X F P = t ^a A/ (i^)n^W M^n^i ( i2 - 4 ) 

ic[p] V iei J \ jeJ J 

where J = [p]\I and, for any function G(t, xi), we denote 

d\ x i\G 

d ^ G = — — (t,l,...,l) 

dx il 1 ■ ■ ■ dx C 

if I = {zi, . . . , i r } with i\ < . . . < i r . Now 

^ f^i(x/)n^) = E(- i ) |Ail " k M/iiip' 

V i€l / cr<Aj iei % ' 

where the sum runs over all non-negative |/| -tuples a = (<Tj)i e / that are less than or 
equal to A/. The second derivative contains fewer terms: 

where the sum runs over all | J|-tuples e = (sj)jeJ such £j G {0, 1} for all j. 

Let us now bravely replace the two derivatives occurring in 1)12 .4|) by their sum- 
expressions given above, and (mentally) expand the product of these sums. This 
gives d\F p as a sum over /, a and e. In this sum, the series d\F p appears twice, 
namely 

(i) for J = 0, a = and e = (0, . . . , 0), 

(ii) for / = \p], a = A and e = 0. 

The corresponding summands are the same in both cases, namely tFod\F p . Hence 1)12.4)1 
can be rewritten as 

(1 - 2tF )d x F p = t SUMMAND, 

I.cr.e 
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where the sum now excludes Cases (i) and (ii). In this sum, all terms of the form 
d T Fk now satisfy k + |r| < p + |A|, so that the induction hypothesis applies to them. 
Note also that (1 — 2t7b) = \f\ — At, so that the previous equation really reads 

VT^tt d x F p = t SUMMAND = RHS. (12.5) 

I,cr,£ 

The latter observation is the key in our proof of (|12.3|) . 

In the right-hand side of the equation, let us study separately the cases where 7 
or J are empty. 

First case: 7 or J is empty. The contribution of the terms for which J = is 

p 

tF ^d X - e F P l[\ £ /. (12.6) 

e^O j=l 

The contribution of the terms for which J = (that is, I = [p]) is 

iF ^(-l)l A l-H 5(Tj p p TJ^:. (12 . 7 ) 

We observe that the terms for which |e| = 1 in (|12.6|) cancel out with the terms for 
which \a\ = |A| — 1 in (|12.7|) . (More generally, the term associated with e, when |e| 
is odd, cancels out with the term associated with a = A — e, but we do not need this 
property). After these cancellations, all the terms <9 T 7fc that appear in this part of 
RHS satisfy k = p and |r| < |A| — 2. In particular, e r < e\ — 1/2 for each of them. 
The induction hypothesis then guarantees that this part of RHS is of the form 



rtriSi — , 

(1 - 4t) e A-V2 

for two Laurent polynomials P\(t) and Qi(t). Given that we still have to divide 
RHS by VI - 4t to obtain the expression of d\F p (see (jl2.5|) ). this part of RHS is 
compatible with the expected form ()12.3|) . 

Before turning our attention to the case ^ 7 ^ [p], let us work out the value 
of Pi (1/4), at least when |A| is even. In RHSi, the only terms d T F p for which 
e T = e\ — 1/2 are those for which |r| = |A| — 2. That is, the terms for which |e| = 2 
in (|T2~oT) . and the terms for which |<r| = |A| - 2 in (|T2~T|) . As 7b -> 2 when t 1/4, 
this means that 



RHSi 



where e^j is the p-tuple that has a one at positions i and j, and zeros elsewhere. 
In the second sum, the partitions a such that C{ = Aj — 2 for some i contribute 
Aj(Aj — 1), while those for which <7j = \ — 1 and aj = \j — 1, so that a = A — Sij, 
contribute AjAj, as in the first sum. A concise way of merging both sums consists in 
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using the notation of (|4.3j) and writing 




RHSi = * I 



(1 - 4t) e A~l/2 
so that the polynomial P\(t) satisfies 

cr<A,H = |A|-2 V 7 

where we recognize the second part of (|4.3[) . 

Second case: 7^0 and J / 0. In that case, the induction hypothesis (|12.3|) 
applies both to S(x-F|/| and d\j- e F\n. Moreover, 



1 

e CT + e Xj ^e <e\ I +e\ J =p+- 



|A/| 



1 

+ 2 



|Aj| 



l<e A -^. (12. 



This implies that the part of RHS for which ^ I ^ \p] can be written as 



dttc P2(t) + Q 2 (t)VT=Tt 
RHSa " (1 - 4t)«*-V2 ' 

for two Laurent polynomials P2(0 and Q2(t)- Given that RHSi has also this form, 
we can conclude at last that 1)12.3(1 holds. 

Let us finally work out the value of P2(l/4), at least when |A| is even. The only 
way for the inequalities (|12,8|) to be equalities is to take a = A/, e = 0, with |A/| 
and |Aj| even. Going back to (|12,4j) . this means that the dominant contribution in 
RHS2 is given by 

4(l-4t) e A-i/2 Yl C ^J- 

In other words, 

which gives the first part in 1|4.3|) . This completes the proof of Proposition 112.11 and 
thus of Theorem 14. II □ 

The proof of Theorem 14.51 is almost the same, using the generating functions and 
recursion relations in Proposition 110.11 and replacing 1 — 4t by 1 — 8t and 1 — 12i, 
respectively. We omit the details. 

To conclude this section, let us sketch the proof of Theorem 14.61 We define the 
generating functions as in (|8.1|) . but replacing £(v) by the depth d(v). Then (|8.2|) 
holds with Xi replaced by Xj. In particular, Fq is still given by (|8.3j) and each Fk(x) 



admits a rational expression in Fq and the Xi. We claim that then (|12.3j) holds, with 
ex =P+ |A|/2- 1/2, and 

dx + o(l) 



dxF p 



(1 _ 4f)p+|A|/2-l/2- 
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This is proved by induction as above. Note that, after p2.4|> . the d\j term expands 
exactly as the d\j term, without (— l)l A d~M an d thus without cancellation; this 
ultimately explains why the exponents e\ increase faster for the horizontal profile 
than for the vertical. The rest is as above. 

The same applies to plane trees, with defined as in (|1U.1|) with £(v ) replaced 
by d(v), and the recursion relation 

F k ( X ) = i+tj2 (n^) ^1(^)^1 ji^j). 

ic[k] / 
where J = [k] \ 1. We omit the details. 
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